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Abstract. We study the shape stabiUty of disks moving in an external Laplacian field in two 
dimensions. The problem is motivated by the motion of ionization fronts in streamer-type electric 
breakdown. It is mathematically equivalent to the motion of a small bubble in a Hele-Shaw cell with 
a regularization of kinetic undercooling type, namely a mixed Dirichlet-Neumann boundary condition 
for the Laplacian field on the moving boundary. Using conformal mapping techniques, linear stability 
analysis of the uniformly translating disk is recast into a single PDE which is exactly solvable for 
certain values of the regularization parameter. We concentrate on the physically most interesting 
exactly solvable and non-trivial case. We show that the circular solutions are linearly stable against 
smooth initial perturbations. In the transformation of the PDE to its normal hyperbolic form, 
a semigroup of automorphisms of the unit disk plays a central role. It mediates the convection of 
perturbations to the back of the circle where they decay. Exponential convergence to the unperturbed 
circle occurs along a unique slow manifold as time < — > oo. Smooth temporal eigenfunctions cannot 
be constructed, but excluding the far back part of the circle, a discrete set of eigenfunctions does 
span the function space of perturbations. We believe that the observed behaviour of a convectively 
stabilized circle for a certain value of the regularization parameter is generic for other shapes and 
parameter values. Our analytical results are illustrated by figures of some typical solutions. 
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1. Introduction. 

1.1. Problem formulation in physical and mathematical context. The 

mathematical model considered in this paper is motivated by the physics of electric 
breakdown of simple gases like nitrogen or argon [H El El IH [5] . During the initial 
'streamer' phase of spark formation, a weakly ionized region extends in a strong 
externally applied electric field. As the ionized cloud is electrically conducting, it 
screens the electric field from its interior by forming a thin surface charge layer. 
This charged layer moves by electron drift within the local electric field and creates 
additional ionization, i.e., additional electron-ion-pairs, by collisions of fast electrons 
with neutral molecules. We here approximate the ionized and hence conducting bulk 
of the streamer as equipotential. In the non-ionized and hence electrically neutral 
region outside the streamer, the electric field obeys the Laplace equation. The thin 
surface charge layer can be approximated as an interface which moves according to 
the electric field extrapolated from the neutral region onto the interface. We therefore 
are concerned with a typical moving boundary problem. 

Such moving boundary problems occur in various branches of physics, chemistry 
or biology. The most extensively studied examples are viscous fingering observed in 
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two- fluid flows ^ or the Stefan problem of solidification from an undercooled melt [7] . 
Other physical phenomena like the motion of voids in current carrying metal films [5] 
lead to similar mathematical models [5] . 

We here discuss the streamer model in two spatial dimensions, where in the 
simplest 'unregularized' version the basic equations coincide with those describing the 
motion of a small bubble in a liquid streaming through a Hele-Shaw cell [TUIITTIIT^ITB] . 
which is a special case of two fluid flow. The unregularized streamer model has 
been discussed in Ref. [H [14] . Restriction to two dimensions in space allows us to 
use standard conformal mapping techniques [6l [15] to reduce the moving boundary 
problem to the analysis of the time dependence of the conformal map that maps the 
unit disk to the exterior of the streamer. 

It is well known that unregularized moving boundary problems of this type are 
mathematically ill posed [15 , in the sense that the moving interface generically de- 
velops cusps within finite time which leads to a breakdown of the model. To suppress 
such unphysical behavior, the models are regularized by imposing nontrivial boundary 
conditions on the interface. For viscous fingering typically some curvature correction 
to the interfacial energy is considered. For the streamer problem a mixed Dirichlct- 
Neumann boundary condition can be derived jTHHGj by analyzing the variation of the 
electric potential across the screening layer. Such a boundary condition is well known 
from the Stefan problem, where it is termed 'kinetic undercooling'. It rarely has been 
considered for Hele-Shaw type problems. There are strong hints [TH [TT] [ISl [IH] but 
no clear proof that it suppresses cusp formation. In particular, it has been shown that 
an initially smooth interface stays smooth for some finite time interval. 

Here we consider the linear stability of uniformly translating circles in a Laplacian 
potential (p that approaches a constant slope (p oc x far from the circle; this means 
that the electric field E = —Vp is constant far from the circle. Though this field 
breaks radial symmetry, uniformly translating circles are exact solutions of the regu- 
larized problem [14j . But perturbations of these circles do not simply grow or decay 
locally as on a planar front or on circles in a radially symmetric force field [T71 [T5] , 
but they are also convcctcd along the boundary; this convection turns out to be a 
determining part of the dynamics. Though physical streamers are elongated objects 
frequently connected to an electrode, the front part of a streamer is well approximated 
by a circular shape. Since it is this part that determines the dynamics, our analysis 
should be relevant also for more realistic shapes like fingers where no closed analytical 
solutions of the regularized uniformly translating shape are known [TH] . In the sequel 
we will use the term 'streamer' to denote the translating circles, being aware that this 
is a slight abuse of the term. 

1.2. Overview over content and structure of the paper. Regularization 
of the streamer model introduces some parameter e that measures the effective width 
of the interface relative to the typical size of the ionized region. The regularized 
problem allows for a class of solutions of the form of uniformly translating circles, and 
linear stability analysis of these solutions can be reduced to solving a single partial 
differential equation. For the special case e = 1, the general solution of this PDE can 
be found analytically, as we briefly discussed in [14]. The present paper is restricted 
to this special case as well. 

The main results of the letter [T^ are the following: The dynamics of infinitesimal 
perturbations are governed by a subgroup of the automorphisms of the unit disk. 
Generically, these automorphisms convect the perturbations to the back of the moving 
body. Initially, perturbations might grow, but asymptotically for time t oo, they 
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decay exponentially. Furthermore, this final convergence back to the unperturbed 
circle follows some universal slow manifold. 

The present paper contains a detailed derivation, discussion and extension of the 
results presented in the letter [14] . Furthermore, the analyticity and completeness of 
temporal eigenfunctions and the Fourier decomposition of perturbations are discussed, 
limit cases of the dynamics are worked out analytically and results are demonstrated 
in a set of figures. 

In detail, the time evolution determined by a PDF is often analyzed in terms of 
temporal eigenfunctions. For the present problem in a space of functions representing 
smooth initial perturbations of the moving circle, no such eigenfunctions exist. They 
can only be constructed if we allow for singularities on the boundary. We here find that 
a subset of these functions with time dependence e""'^, n G Nq, is intimately related to 
the asymptotic convergence of the perturbations. These functions show singularities 
only at the backside of the circle, and the front part of any smooth perturbation 
can be expanded in this set of functions. The spatial domain of convergence of this 
expansion increases with time and asymptotically for t — > oo, it covers almost the 
whole streamer. In this restricted sense these eigenfunctions form a complete set. 

These results dealing with infinitesimal perturbations, of course, do not imply the 
asymptotic stability of the circular shape against finite perturbations. To solve this 
problem, the full nonlinear theory must be considered. Nevertheless, a first hint might 
be gained by considering the evolution of a finite perturbation under the linearized 
dynamics. Due to the conformal mapping involved, the absence of cusps under this 
evolution is not a completely trivial question. We here show that for a large range of 
smooth initial conditions, the shape of the streamer stays smooth under the linearized 
dynamics. 

All the present work deals with the exactly solvable case e = 1 whereas the phys- 
ically most interesting case is e <C 1. We, however, believe that the features we could 
identify explicitly for e = 1, are generic for all e > 0. In particular, the subgroup of 
automorphisms of the unit circle leads to the basic mechanism of convective stabi- 
lization, it is for all e > intimately related to the characteristic curves of the PDE, 
and it also governs the dynamics in another exactly solvable case, namely for e = oo. 
Furthermore, it can be shown [W that the temporal eigenvalues A„(e) emerging from 
A„(l) = — n, stay negative for all e > 0, which also indicates that the circle might be 
asymptotically stable for arbitrary e > 0. 

This paper is organized as follows. In Sect. [5] we introduce the model, and the 
linear stability analysis of translating circles is carried through in Sect. [31 These two 
sections are extended versions of Ref. [Hj. Analytical results based on the PDE of 
linear stability analysis are derived in Sect. [H in particular, center of mass motion, 
internal motion, (non)analyticity and completeness of eigenfunctions, intermediate 
growth and asymptotic decay of perturbations, Fourier representation and motion 
of nonanalytical points in the complex plane of the conformal map. These dynamic 
features are illustrated by explicit examples in Sect. [5] The appendix contains a 
discussion of the case e = oo. 

2. Physical model and conformal mapping approach. 

2.1. The model. We assume the ionized bulk of the streamer to be a compact, 
simply connected domain Vi of the (x, ?;)-plane. Outside the streamer, i.e. in the 
open domain T^n, there are no charges and the electric potential obeys the Laplace 
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equation 

(2.1) Aif^O for e P„. 

The streamer moves in an external electric field that becomes homogeneous far from 
the ionized body; therefore the electric potential ip at infinity obeys the boundary 
condition 

(2.2) If EqX + const for \/ ^ oo . 

This condition excludes a contribution to 95 diverging as ln(a;^ + y^) which implies 
that the total charge due to the sum of all electrons and ions vanishes within and 
that the far field has the form 

E = -Vip -Eox, 

where x is the unit vector in a;-direction. On the surface of the streamer we impose 
the boundary condition 

(2.3) ^ = £n-V^, 

where fi is the unit vector normal to the surface pointing into 2?„. Here as well as in 
Eq. (12. 4p below it is understood that the surface is approached from T>n ■ As mentioned 
in the introduction, this boundary condition results from the analysis of the variation 
of the potential across the interface, and the length parameter £ can be interpreted 
as the effective thickness of the screening layer. The case i = corresponds to the 
unregularized case with a pure Dirichlet condition on the moving boundary. Dynamics 
is introduced via the relation 

(2.4) v„=n-yip, 

which holds on the boundary and determines its normal velocity Vn- This defines our 
model. For further discussion of its physical background, we refer to [Tl [2l [3l l4l [5l fT6] . 

Now obviously, Eq can be absorbed into a rescaling of the potential (p and of the 
time scale inherent in the velocity ?;„, therefore henceforth we take Eq = 1. Clearly 
the model defined here is most similar to a model of the motion of a small bubble in 
a Hele-Shaw cell [UlIIl], except that the boundary condition (|2.3p is of the form of 
a kinetic undercooling condition [T71 [TH] . 
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2.2. Conformal mapping. A standard approach to such moving boundary 
problems proceeds by conformal mapping [HI HI]- We identify the (x,?/)-plane with 
the closed complex plane z = x+iy, and we define a conformal map f(uj,t) that maps 
the unit disk Uui in the w-plane to I?„ in the z-plane, with w = being mapped on 

z = oo 



(2.5) 



+ f{LO,t), a_i(0>0. 



Here the function / is holomorphic for uj G U^, and we assume that the deriva- 
tives d" of all orders n exist on the unit circle dUu,- This restricts our analysis 
to smooth boundaries of the streamer. (Weaker assumptions on boundary behavior 
briefly will be discussed in Sect. 14. 8h We recall that the closed physical boundary can 
now be retrieved as Xa{t) — 3?/(e*",i) and ydt) — 3/(e*",i)) where the interface 
parametrization with the real variable a £ [0, 2t:[ is fixed by the conformal map. 

By virtue of Eq. (j2.ip . the potential ip restricted to is a harmonic function, 
therefore it is the real part of some analytic function $(z, t), which under the conformal 
map (j2.5p transforms into 



(2.6) 



a-i{t) 



+ $(w,i)- 



Here the holomorphic function $ obeys the same conditions as / above, 
results from the boundary condition (|2.2p with Eq — 1, and (|2.5p . 
Conditions and ((^ take the form 



The pole 



(2.7) 



dtf 



for UJ e dUu: 



for UJ £ dUrj 



Eqs. (|2.5p - (|2.8p form the starting point of our analysis. 

3. Linear stability analysis of translating circles. 

3.1. Uniformly translating circles. A simple solution of Eqs. (|2.7p 
takes the form 



f(^){uj,t) 



(3.1) 



R 


2R 


t, 


UJ 


^ R + £ 




'I R 


-I 


R 




f 


UJ R- 



In physical coordinates x and y, it describes circles of radius i? > centered at 
x{t) — vot and moving with velocity vq — 2R/{R + £) in direction x. Thus the point 
UJ — 1 maps to a point at the front and the point uj ~ —1 maps to a point at the back 
of the streamer. These points will play a crucial role in our analysis. 

We note that the one-parameter family (13. ip of solutions parametrized by R, that 
is found in the regularized model, is a subset of the two-parameter family found in 
the unregularized case i = 0. As is well known, for £ = all ellipses with one axis 
parallel to x are uniformly translating solutions [TUj . 
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3.2. Derivation of the operator for linear stability analysis. We now 

derive the equation governing the evolution of infinitesimal perturbations of the circles 
(j3.ip . In general, the parameter R can become time dependent. We use the ansatz 



(3.2) 



/(^,0 



R(t) 



- x{t) +r] (3{Lu,t) 



< <P{uj,t) = R{t) 
dtx{t) 



1 R{t) - e 
u ^ R{t) + a 



cj + 77X(cj,t) 



R{t) + 1 



where (3 and x ^re holomorphic functions of w and where ry is a small parameter. 
However, working to first order in rj it is found that R stays constant. This results 
from the fact that the dynamics embodied in Eq. (|2.8p strictly conserves the area \T>i\ 
of the streamer, which in this context is the equivalent to the temporal conservation 
of the zero order Richardson moment [T31 [HI [H] , but integrated over the complement 
of P„. In terms of the mapping /, the conserved area {DA can be written as 



|A| = da{^[f{e'^,t)\ -a;(t))5„5[/(e'",i)] 



(3.3) 



7ri?2(t) -^jf j da^ [/3(e*", t)] 9^,3 [/3(e*", i)] 



Now introducing the time independent length i?o through |2?i| = 7ri?g, we find R{t) = 
i?o+C(??^), which proves that R is time independent within linear perturbation theory. 
In the sequel we will use i?o as our length scale, introducing 



(3.4) 



e — — and 
Ro 



t 



1 + eRa' 



and rescaling / and $ by factors I/Rq. We note that within a dimensionless time 
interval r of order unity, the streamer moves a distance of the order of its size. 

With the thus simplified ansatz (|3.2p . Eqs. (|2.7p and (|2.8p evaluated to first order 
in 7] take the form 



(3.5) 



3? 



Lu{du, - dr)f3 - ^-^^ujdu,x 



= 0, 



for uj G dUu; 



)R [e(L^2 ^ i)^Q^p - (1 + e)(l + euod^)x] = , 
Since [3 and x are holomorphic for lj G lA^^, these equations imply 

(3.6) 



for UJ eU^ 



e{uo^ + \)ojd^l3 - (1 + e)(l + eojd^)x = ia{t) , 



where a{t) is some real function of time, x is eliminated by substituting the expressions 
for di^x and d^x from the first equation and its derivative into the second equation 
differentiated with respect to uj. This yields 



(3.7) 
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where is the operator 

(3.8) = I 9^ (c^' - 1)^ + e ujd^dr + (1 + e) 9, - 9^ . 

3.3. Normal form of and induced automorphisms of the unit disk. It 

is instructive to transform to the normal form of a hyperboHc differential operator. 
We introduce 

(3.9) r = tanh^, 

mapping the time interval r S [0, oo[ to T G [0, 1[, and 

cj + T 



(3.10) C = 

to find 



(3.11) = th{C,T)dTd^ + + (1 + e)^T , 



9MC,r) 

(C-T)(1-TC) 



(3.12) 1-T^ • 

This identifies the manifolds T — const or Q — const as the characteristic manifolds 
of our problem for all e 7^ 0. 

As function of the 'time- like' parameter T, < T < 1, the transformation = 
C(a;,T) in Eq. p.lOp represents a semigroup of automorphisms of the unit disk, with 
fixed points 

Q = UJ = ±\. 

For T— >1, corresponding to r^oo, all points ijj ^ —1 are mapped into C = +1, so 
that the large time behavior of any perturbation is governed by this attractive fixed 
point. 

3.4. Analytical solutions of Eq. (13. 7p for special values of e. The general 
solution of Eq. (|3.7p can be found analytically for the special values e = 0, e = ±1 
and e = cxD. In the unregularized case e = 0, evidently any function 

/3(c^,r)=/3(c^ + r) 

is a solution, and any singularity of /? found in the strip 

< < 00, - 1 < < 1 

will lead to a breakdown of perturbation theory within finite time. This is the finger- 
print of the ill-posedness of the problem for e = 0. 

For e — —1, /3(w,t) generically for all r > has a logarithmic singularity at 
LL> — —T(t). We recall that negative values of e = £/-Ro imply negative thickness of 
the screening layer and thus are of no physical interest. 

The case e = +1 is discussed in detail in the remainder of the paper. Though 
a regularization length t identical to the object size i?o is somewhat artificial, it is 
accessible to rigorous analytical treatment and, as explained in Section 1.2, we expect 
it to reveal generic features of the behavior for all e > 0. 

This is supported by the results for e = 00 which show essentially the same 
features as the results for e = 1 below. Though the limit e — > 00 is physically absurd 
when applied to streamers, it is worth studying with respect to the properties of the 
operator >Cc, and we present a short discussion in the appendix. 
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4. Strong screening: analytical results for e — 1. 

4.1. Analytical solution of the general initial value problem. With the 
form ((3111) of the PDE dSJ]) for e = 1 reduces to 

(4.1) dT [2 + hiC,T)dc) P^O, 
showing that the function 

(4.2) G(C)-(2 + MC,m)/3 

is independent of T. To determine (3, we use Eq. (|3.12p : h{(,T) = uj/d^uj to find 

(4.3) (2 + c.a^)/?(c.,T) = G(C), C = C(c^,T(t)). 
The solution regular at w = takes the form 

(4.4) /3(.,r). + 



1 + xT{t) 
u 

A second independent solution is singular in w = 0: 

(4.5) /3,ing(cj,r) = ^. 

The function G in the regular solution (j4.4p is determined by the initial condition 
/3(w, 0) through 

(4.6) G{uj) = {2 + ujdu;) f3{uj,0). 

It thus is holomorphic for w in the unit disk and all derivatives exist on dLi^j, since 
we assume the initial surface to be smooth. Eq. (|4.4[) then shows that /3(w, r) inherits 
these properties for all t < oo. 

4.2. Automorphism of unit disk and a bound on the perturbation. It 

is now clear that the automorphisms ^(w, T) of U^i from Eq. (|3.12p contain the basic 
dynamics and, as shown in the appendix, this also holds for e = oo. This is to be con- 
trasted to the unregularized case e = 0, where the dynamics amounts to a translation 
of the unit disk. With the present dynamics, in the course of time larger and larger 
parts U{S) of the unit disk are mapped to an arbitrarily small neighbourhood 
IC — 1| < (5 of the attractive fixed point C = 1. According to Eqs. (|4.4p and (14. 6p . 
the initial condition in the neighbourhood ja-" — 1| < S then determines the evolution 
of (3{uj,t) in all U{6). As a consequence, any pronounced structure found initially 
near ojq, \ujq — 1| > (5, is convected towards w = — 1. Quantitatively this behavior is 
embodied in Eq. (|4.17p below, and explicit examples will be presented in section [SJ 
see, in particular, figure !?^ . 

For the further discussion we normalize G(w) so that 

(4.7) max |G(cj)| = 1. 
Eqs. (|13), ((i?71) yield a boimd on I3(uj,t): 

(4.8) |/3(cc>,t)| < i; |c^| < 1 , 0<T<1. 

Thus the perturbation can shift the position of the streamer at most by rj/2, and 
therefore it cannot affect the asymptotic velocity of the propagation. 
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4.3. Center of mass motion for < r < oo. In precise terms tlie position of 
the streamer can be defined as tlie center of mass 

(4.9) Zcm = Xcm + Wcm = J dx dy {x + iy) , 

wliere the integral is related to the first order Richardson moment. Evaluating 
Eqs. (|4.9p and (|4.4p . we find to first order in 77 

(4.10) Zcm = r + ??/5(0,T), 

Here r is the uniform translation of the unperturbed circle. The additional center 
of mass motion (|4.1ip for all times is explicitly given by the initial condition /3{uj,0) 
through Eq. (14.611 and the transformed time variable T{t) from Eq. (j3.9|l : for t 00, 
it approaches /3(0, r) G(l)/2. 

4.4. Internal motion: convergence along a universal slow manifold for 

T — > 00. We now concentrate on the perturbation of the circular shape, given by 

(4.12) PiLO,T)=PiLO,T)-p{0,T). 

The explicit expression 



(4.13) /?(c^,r) 







GiT) 



yields 

(4.14) lim (3{uj,t) = 



for arbitrary G, i.e. for arbitrary initial condition (j4.6|) . Thus the shape perturbation 
converges to zero as r — > 00, and the circular shape is linearly stable. 

We note that this holds despite the fact that the limits ui ^ —1 and t 00, (i.e. 
T ^ 1), do not commute 

lim lim G(C(w,T)) = G(-l), 
lim lim G(C(c^,T)) = G(+l). 

This peculiar behavior near the backside of the streamer, at co = —1, shows up only 
in the rate of convergence. 

Investigating the rate of convergence for r — > cxd, we first exclude a neighborhood 
of a; = — 1 and expand G in the integral (|4.13|1 as 

where G'(w) = 5^G(w) . 

With 
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the integral yields 

/3(a;,T) 



(4.15) 
valid for 



G"(l) 



ln(l 



ijj 



+ 0(e-'") , 



1 + w| > Iwle" 



Thus outside the immediate neighborhood of w = —1, the shape for all smooth initial 
conditions with G"(l) ^ convergences exponentially in time as e""^ along a universal 
path in function space, given in Eq. (|4.15[) . For G"(l) = the first non- vanishing term 
in the expansion of G dominates the convergence. 

To analyze the neighbourhood of w = — 1 we take the limit r ^ oo, with 



(4.16) 

fixed. We find 

/3(C^,T 



s = (1 + cj)e^ 



G'(l) 



= 4 (ln(2 + s) - t) e-^ 



2G'(1) +41n 



G' 



s-2 



G 



-G'(l) 

4/(2+s) 



G'(l) 



+ (2 + s) (^G(l 
(4.17) +0 . 

In terms of tj, the first contribution on the r.h.s. takes the form 
4 (ln(2 + s) - t) e"^ = 4e"^ In (2e"^ + 1 + 



which shows that a logarithmic cut of /3(w, t) reaches u; = — 1 for t — > cxd, but with a 
prefactor vanishing exponentially in that limit. We thus have found a week anomaly 
of the asymptotic relaxation near w = — 1: In a spatial neighborhood of order e^"^ the 
exponential relaxation is modified by a factor r. Furthermore, as mentioned above, 
all the initial structure of /3(a;, 0) is compressed into that region. This is obvious from 

the occurence of G (f^) etc., in Eq. (|4.17p . 

To summarize, we have found that the shape of the interface for r ^ oo converges 
to the circle along a universal slow manifold (I4.15p . except for a weak anomaly (|4.17p 
at the backside at a; = — 1. 

4.5. (Non-)analyticity of temporal eigenfunctions. In many cases, a full 
dynamical solution for arbitrary initial values cannot be found, and rather temporal 
eigenfunctions are searched for. However, in the present problem, functions f3{L0,T) 
resulting from smooth initial conditions cannot exhibit exponential behavior in time 
for all r, < T < oo. This is seen easily by introducing 



(4.18) 



G{x) = G 



x-l 
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writing G(C) in the equivialent form 



I+ujTJ \u} + 1 

and substituting this form into Eq. (|4.4p . Postulating strict exponential time behavior 
/3 ~ e^^^ one finds 



(4.20) /3(o.,T)oce^^/3AM, f3x{uj) = I dpp 



ojp — 1 



ujp + I 



Any eigenfunction P\{lu,0) with A 7^ clearly is singular either at w = +1, or at 
a; = —1, or at both points. It therefore conflicts with smooth initial conditions. On 
the other hand, omitting a neighbourhood of — —1, eigenfunctions exist for all 
-A e No. 

4.6. Completeness of the eigenfunctions near w = 1. In some neighborhood 
of w = 1, we even can show that any regular solution (3{uj, t) can be expanded in terms 
of the 'eigenfunctions' /3_„(cl;), n G Nq. This results from the Taylor expansion 

00 

(4.21) G(y) = ^g„y", 

n=0 

which by assumption converges in a disk of radius f > 0: Rewriting Eq. (j4.4p as 

^, , f^xdx^/x + T\ f^xdx^fx + T 
P{'^,t)= -^G{— — -]- G 







1 + xTj /,, ^2 \l + xT 



u^^. M{T) ['xdx (l~x 

and (3-n{u>) in a similar form as 

Mn x dx f \ ~ X 



we find 



R ( \ 



(4.24) /3(^,r) = ^ + f;g„ 

Provided < f, we can separate the sum into the contribution oc l/w^ and the rest. 
Since both (3{ijJ,t) and the /3_„((jj) are regular at w = 0, the contributions oc l/w^ 
have to cancel, which yields the final result 

oc 

(4.25) /3(c.,T) = ^5„/3_„(a;)e-"^ 

Tl = 

This result is valid for e^'^ < f in the disk 
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/3(1,t) 




Fig. 4.1. P{e^°' , r) from Eq. 14.2T^ for a = as a function of subtracted time 9 = t — ln2fc. 



It generalizes the asymptotic result (|4.15l) . Indeed, the universal shape relaxation 
found in (|4.15p together with the center of mass relaxation (|4.1ip precisely follow the 
slowest eigenfunction from (|4.2ip with A — —n = — 1. Furthermore this result shows 
that the range of validity of the expansion (|4.25p increases with r and asymptotically 
covers the whole complex plane except for the special point u) = ~\. 

4.7. Intermediate temporal grovifth and coupling of Fourier modes. 

Having found that the space of regular functions does not allow for strictly exponential 
time behavior, we now consider the typical time variation of smooth perturbations. 
Before the exponential relaxation sets in, such perturbations typically will grow, and 
this growth can be quite dramatic. As an illustration we consider a perturbation 
defined by 



corresponding to initial conditions 



fc > 1 , 



(4.26) 



/5(^,0) 



fc + 2 



For T — 1 — e ®/fc, corresponding to times t = 9 + ln(2fc) + 0{l/k), we can write 



G 



U! - 



1 +wT 



l+uj k 



1 - 



1+w fe 



V 




1 = exp 


1 + 0; 





l + O 



where we again excluded some neighborhood of w = — 1 . Substituting this expression 
into Eq. (|4.13p we find on the unit circle uj ~ e*": 

(4.27) /3(e'",r) 
1 

2ipsina 



dpp exp 



P 



1 



2p cos a 



■ exp [- 



O 



Figure 14711 shows this function, evaluated at a = (tj = 1). The behavior is quite 
peculiar. Up to times of order Infc the perturbation stays of order 1/fc ^ 1, then it 
increases roughly exponentially up to values of order 1, and finally it decreases again 
exponentially, approaching the slow manifold (|4.15p . Thus for very large fc the initial 
perturbation P(uj, 0) ^ 1/fc in some time interval can be amplified by a factor of order 
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fc, and Eq. (|4.27p shows that the leading behavior in that time interval is independent 
of fc. 

Closer analysis shows that in terms of a formal Fourier expansion 

oo 

(4.28) /3(e'",T) = ^a„(r)e™", 

n=l 

the amplification is carried by the low modes, n — 0(1). As will be illustrated by an 
explicit example below, cf. figure \5.2b . in such a mode representation the time evo- 
lution feeds the strength of the perturbation sucessively into lower and lower modes. 
This is equivalent to the observation that the automorphism e*" — * C(e*",r) drives 
all the pcrturbative structure towards a = tt and smoothens the remainder of the 
interface. Note, however, that starting with a perturbation ~ uj'' in the course of 
time also modes n > k are (weakly) populated to build up a complicated structure 
near uj = —1. We recall that for the unrcgularized model e = 0, the time evolution of 
a perturbation oc w'' populates only modes fc < n . 

4.8. Motion of the zeros of d^f and cusps. So far we have shown that the 
propagating circle is linearily stable, i.e., we implicitly considered perturbations of 
infinitesimal strength rj. The full nonlinear evolution of a finite perturbation is beyond 
the scope of this paper. Still, it clearly is a question of practical interest, whether a 
finite perturbation evolving under the linearized dynamics, for all times satisfies the 
assumptions underlying the conformal mapping approach. For the mapping to stay 
conformal, all the zeros of duif{oj,T) must stay outside the unit circle. We thus here 
analyze the roots of the equation 

(4.29) = 5^/(c^,T) = -^+77 9^/?(a;,T). 
Using Eqs. (|4.3p . (14. 4p . we can rewrite this equation as 



1 



(4.30) 2r] J dpp 







1 + LuT J \1 + pujT J 



With our normalization (|4.7p of G, for all uj in the closed unit disk the l.h.s. of this 
equation is bounded by 2\ri\. We conclude that the bound 

(4.31) \v\ < \ 

guarantees that within the framework of first order perturbation theory the mapping 
stays conformal for all times. We now will show that this bound in general cannot be 
improved. 

For r — > oo, zeros of du)f{uj, t) reach uj — —I, which is a consequence of the fact 
that in this limit an infinitesimally small neighborhood of a; = — 1 under the mapping 
w — >• ^ is mapped essentially on the whole complex plane. We now analyze this limit 
for the simple example G(w) = u. Substituting this form into the asymptotic behavior 
(|4.17p and using the definition (|4.16|1 of s, we find 
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Eq. (|4.29p reduces to s = 4r/ — 2, showing that a zero luq of 9^/(a;,r) approaches 
w = — 1 as 

wo = -l + (4??-2)e-^. 
For clIq to come from outside the unit circle we clearly must have 

(4.32) 5RM < 1 . 

To get some feeling for the estimate (|4.3ip . we note that for G{uj) = uj*' the map 
initially (for r = 0) is conformal provided |?7| < 1 + 2/fc. We conclude that under the 
linearized dynamics a large part of smooth initial conditions relaxes to the circle. 

In all the discussion of this section we have assumed the initial boundary to be 
smooth, so that all derivatives d^G{io) exist on the boundary \w\ — 1. Inspecting the 
results it is obvious that this assumption can be considerably relaxed, since only those 
derivatives which show up explicitly, have to exist. Thus, for exponential relaxation 
(|4.15p outside the neigborhood of = — 1 to prevail, the existence of dujG{uj) is 
sufhcient, which amounts to the condition that the curvature of the initial boundary 
is well defined. For the circle to be linearly stable, Eq. (|4.14p . it is sufficient that G(e*") 
is bounded and continuous, which implies that the boundary has a well defined slope. 

If the initial boundary shows a cusp, the time evolution sensitively depends on 
the details. If the cusp is found in forward direction, so that G{uj) diverges for 
u! ^ 1, the streamer will be strongly accelerated. In a related model [12], such an 
effect has been pointed out before. Furthermore the shape will not relax to a circle, 
and the conformal map will presumably break down at finite time. If the cusp does 
not affect the analyticity of G(u;) near w = 1, it is convected towards the back and 
broadened, whereas the front of the streamer approaches the circular shape. Still, 
however, conformality of the map may break down at finite time. 

5. Explicit examples for e = 1. We here illustrate the general results by some 
examples. 

5.1. The evolution of Fourier perturbations. We first consider perturba- 
tions of the form 

(5.1) /3W(^,o) = ^, i.e.G{Lo)=co^ 

The integral (|4.4p is easily evaluated to yield 

P^'Hu^.r) = _i-^|r^-+((Ta;)2_i) ^'^ 



-fc(i-r 



2\ 



r^--(.r+i)c^- + l+^+^'(i-^) 



(5.2) -(Ml+cT)-!]— (r-T^ 



^ 1 rpi^ 



V 

v=l 



where T = T{t) and C, — Ci'^, T{t)) are given by Eqs. (|3.9p . or (|3.10p . respectively. 
In figure 15.11 we have plotted snapshots of the resulting motion of the interface, de- 
termined as 



(5.3) z ^ X + iy ^ - +T + ri f3^''^LLi,T) , ui 



e'" , < a < 27r . 
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Fig. 5.1. Snapshots of the evolution of the streamer for k = 2 (left column) and fc = 10 (right 
column) at the indicated instants of time. The solid lines represent the perturbed interfaces. The 
gray disks move with the center of mass velocity of the perturbed circles. One gray disk has 

been omitted for clarity. See the text for further discussion. 



The direction of motion, i.e. the positive x-direction, is downwards. Together with 
the moving interface we show the unperturbed circular streamer at different times as 
gray disks with the center moving according to 

(5.4) zUr) = r + I G{T{t)) = r + | tanh'^ ^ , 

as predicted for the center of mass motion for the perturbed streamer in Eq. (j4.10p . 

In figure [nH] we perturbed the circle by ry fc = 2orfc = 10, using the same 
parameter = 0.6e™^^ in both cases. The starting position for A: = 10 is shifted 
relative to that for fc = 2 by a distance corresponding to Ar = ln5. As discussed 
below Eq. (|4.27p . for 1 < fci < fc2 we expect 

/3['^il [u, t) » /Jl^'-^l {lu, t + ln(fc2/fci)). 

Figure ISTTl illustrates that such a 'universality' for the gross structure holds down to 
very small fc. (Of course the choice of differing values of rj would distort the figures 
and mask this feature.) Basically during time evolution the initial maximum closest to 
the forward direction is smeared out and builds up the asymptotic circle, whereas all 
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X -1.1 -1.05 -1 5) 0-4 0.6 0.8 rp 

Fig. 5.2. a) Magnified plot of the backside of the streamer for k = 10, rj = 0.6 e"^/'' (as in the 
right column of fiaure \5. l\l for the t values given. The overall motion is subtracted. We observe the 
compression of the fine structure and the intermediate growth of the perturbation. Asymptotically 
for T —* oo, the structure converges to the gray circle. In the comoving frame, the gray dot marks 
X + iy = —1 which is the point to which the structure finally is contracted. Note that the scale of x 
is stretched compared to that of y, and that the figure is turned relative to figure [5J\ 
b) The amplitudes an as in Eg. as a function of T for k = 10; the values of n are given. 




Fig. 5.3. Motion of the zeros of dujf in the u-plane for k = 2 and rj = 0.6 e*'^/'* (as in the left 
column of fiaure \5. The dots give the position for r = 0, 1, 2. The horizontal line is the cut for 
T = 2.51, where one zero enters the second sheet (broken curve). The unit disk is also shown. 

other structures are compressed at the backside. For fc = 10 the comphcated structure 
at the back is magnified in figure [5T2b . Figure [5T2b shows the time dependence of the 
coefficients a„ of the low modes e™" in the expansion (|4.28[) . again for k = 10. It 
illustrates how the strength of the perturbation cascades downwards in n and increases 
in time, until it completely is absorbed into the lowest mode, i.e., the overall shift of 
the circle. We should recall, however, that also modes n > k are weakly populated to 
build up the structure at the back. 

For fc = 2, figure [?751 shows the motion of the zeros of 9^/(a;, r) in the complex lo- 
plane, as discussed in section It corresponds to the fc = 2 part of figure [5TT1 Two 
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zeros, which initiaUy are close to the backside of the unit circle, for t — s- cx) approach 
a; = — 1. They clearly are associated with the two maxima that in the comoving 
frame are convected towards z = x + iy = —1. The third zero, originally found close 
to w = +1, after a large excursion leaves the physical sheet at time r ~ 2.51. The 
logarithmic cut is on the negative axes, with the branchpoint cjbp = ^^/T(t) reaching 
to = —1 for T oo. 

5.2. The evolution of localized perturbations. We finally consider some 
more localized perturbation, defined by 



(5.5) 



(1 - 7)e'"° 
G(^) = — — —r > 7 > 1 , 



uj — je" 

corresponding to an initial perturbation 



(5.6) 



rj I3{bj, 0) = 1] 



(1 _ ry)j 



_,2mo 



In ( 1 - - e"'"« 
7 



— e 

7 



The result for (3{uj,t) reads 



(5.7) P{uj,t) 

where 

(5.8) 



(7 - l)e'' 



r(r) 



^e-iao _ t{t) {2b{T) 



b{r) = 



T{t)\ ln(l + 5(r)w) - fe(r)tj 



b{r) 



1 - T(r)7e^"" 
T(r) - 7e*"o 



We note that ^(t) 1 for T{t) 1, so that in the large time limit the logarithmic 
cut reaches u — ~1. As discussed in the context of Eq. (|4.17|1 . this is a generic feature 
of the present problem. Our choice of parameters 7 — 1.1, ao = — 7r/12, rj — 1.5, 
almost produces a cusp in the initial condition: the only zero of dujf{uj, 0) is found 
at ujQ = 1.001 exp(—.243i). This zero, however, is driven away from the unit circle 
and leaves the physical sheet. Another zero that entered the physical sheet somewhat 
earlier, for t 00 reaches u = —1. Figure [5^4h shows snapshots of the time evolution 
of the perturbed interface in a representation like figure 15.11 It illustrates how the 
peak rapidly is smeared out and the interface becomes smooth. Figure lS^ follows the 
evolution of the peak for short times and shows how it is convected and broadened. 

We finally note that in the special case, where the initial peak strictly points 
in forward direction (ao = 0), convection cannot take place. The peak simply is 
broadened and vanishes, whereas some new peak shows up at the back for intermediate 
times. 

We acknowledge helpful and motivating discussions with F. Brau, A. Doelman, 
J. Hulshof, H. Levine, L.P. Kadanoff, S. Tan veer and S. Thomae. 

Appendix A. The limit e 00. 

For e -> 00, the PDF ([SJ]) with the form ([3?TT|) of jC^ reduces to 
(A.l) (mC,T)<9c + i) aT/3(C,r)-0, where /3(C,T)=/3(c.,t). 
Fq. (jA.ip allows for a large set of solutions obeying the same initial condition 



(A.2) 



/3(l^,0)=/3o(w), 
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but imposing regularity on the unit disk lAui , we single out the simple form 
(A.3) /3(c.,r) =/3o(C). 

Thus for e = oo, the dynamics is simply given by the automorphisms ui > C,{ijJ,T). 

This implies that I3{uj^t) is bounded uniformly in t as 

(A.4) |/3(t^,r) I < max | /3o(^) |, 

so that in contrast to the case e = 1, there is no intermediate growth of the pertur- 
bations. 

The shift of the center of mass is given by (cf. Eq. (|4.10p ): 

(A.5) /3(0, r) = /3o(r(T)) = /3o(l) - 2 /3^(1) e"- + O (g-^^) , 

and except for the point lo — ~\, the shape again converges exponentially in time to 
the circle along the universal slow manifold 

(A.6) t) - /3(0, r) - /3^(1) -i^ e'^ + O {er^^) , 

cf. Eq. (|4.15p for e — 1. Again the neighbourhood of a; = 1 for time r = 0, more 
precisely /3o(l) Sind Pq{1), determine the long time convergence. Since by assumption 
(3o{uj) is analytical at a; = 1, evidently an eigenfunction expansion in the sense of 
subsection 14.51 exists. 

The only major difference to the case e = 1 concerns the point lu = —1. Clearly, 

(A.7) /3(-l,r)EE/3o(-l) 

independently of r, and indeed for t ^ oo the conformality of the mapping breaks 
down in the neighbourhood of a; = —1 since 9i^/3(a;, ■'")|;^__]^ diverges. 
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